Amine-Gold Linked Single-Molecule Junctions: Experiment and Theory 
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The measured conductance distribution for single molecule benzenediamine-gold junctions, based 
on 59,000 individual conductance traces recorded while breaking a gold point contact in solution, 
has a clear peak at 0.0064 Go with a width of ± 40%. Conductance calculations based on density 
functional theory (DFT) for 15 distinct junction geometries show a similar spread. Differences in 
local structure have a limited influence on conductance because the amine-Au bonding motif is 
well-defined and flexible. The average calculated conductance (0.046 Go) is seven times larger than 
experiment, suggesting the importance of many-electron corrections beyond DFT. 

PACS numbers: 71.15.-m, 72.10.-d, 81.07.-b, 82.45.Yz 



Discovering the anatomy of single-molecule electronic 
circuits, in order to exploit their transport behavior, 
poses fundamental challenges to nanoscience Single- 
molecule manipulation requires probes at the limits of 
experimental resolution [J 0, 0, M, H, 0| , as control at 
the single-bond scale is necessary for reliable transport 
properties. The success of amine-Au links for the real- 
ization of single-molecule junctions has allowed system- 
atic measurements of the relationship between molecular- 
scale structure and conductance for several families of 
molecules [1, S, 0| • Experimental trends were explained 
based on the hypothesis of selective bonding between the 
amine lone pair and an undercoordinated gold atom Q. 
Understanding the essential junction features that lead 
to reproducible electrical properties would enable gener- 
alization to other metal-molecule link chemistries, while 
also furthering the development of a predictive theoreti- 
cal approach to nanoscale transport. 

Most recent theoretical investigations of nanoscale 
transport rely on a Landauer approach, simplified to 
treat electronic interactions at a mean-field level within 
density functional theory (DFT). While this framework 
has proven relatively accurate for metallic point con- 
tacts [10|, the computed conductance often substan- 
tially exceeds measured values for molecular junctions 
@- EH, Ell, 13, [3], raising questions about fundamental 
theoretical approximations. However, since the predicted 
conductance can be very sensitive to contact geometry 
(changing by orders of magnitude for thiol- Au junctions, 
for instance T3, 13, lH), uncertainty over junction struc- 
ture has obscured this issue. A clear benchmark is re- 
quired. 

In this Letter, we use the reproducible conductance of 
amine-gold linked single-molecule circuits [1, H, @] to de- 
velop a detailed, atomic-scale understanding of a proto- 



typical molecular junction. A large experimental data set 
for benzenediamine (BDA) is analyzed in depth, and di- 
rectly compared with DFT-based transport calculations. 
We show that amine-gold bonding is highly selective, but 
flexible, resulting in a conductance insensitive to junc- 
tion structure. The computed conductance has a narrow 
spread consistent with experiment, although the aver- 
age calculated conductance is ^7x larger. Our bench- 
mark comparison demonstrates the importance of many- 
electron corrections beyond DFT. 

To obtain a detailed picture of junction structure that 
can be connected to its transport properties, we ana- 
lyze the distribution of low-bias conductance of Au-BDA- 
Au junctions measured by repeatedly breaking Au point 
contacts in a molecular solution [6j. Many conductance 
traces reveal a step below 1 Go (2e 2 /h) (Fig. 1(a)), 
corresponding to a single BDA molecule between the 
two atomic point-contacts 0, Fig. 1(b) shows the 
histogram (black) computed from all 59,000 measured 
traces, from 4 different tip/sample pairs, without any 
data selection. The most probable junction conductance 
obtained from Lorentzian fits to the histogram peak 
is 0.0064 ± 0.0005 G , with half width 47 ± 8 %. 

Not all measured conductance traces have a step at 
a value near the molecular peak in the full histogram 
(black, Fig. 1(b)). Applying an automated sorting algo- 
rithm to the full data set indicates that a molecule is 
trapped between the breaking Au contacts in about 36% 
of the traces. Our algorithm also determines the length 
and slopes of the molecular steps, whose distributions are 
shown in Fig. l(c-d). We find that the majority of steps 
(~75%) are < 1 A long, with a power law distribution 
extending to about 4 A. Further, ^75% of the steps are 
flat, with a slope ranging from -0.15 to -0.05 Go/nm. 

Further insight is gained from comparing the distri- 
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FIG. 1: (color online) (a) Sample conductance traces. 
Dashed lines are linear fits to the molecular step determined 
by an automated algorithm (applied bias: 25 mV). (b) Con- 
ductance histograms of all 59000 measured traces (black) and 
of 21312 selected traces with steps (red), shown on a log-log 
scale. Inset (top-right): Conductance histogram of selected 
traces (red) and histogram of average step conductance (blue) . 
Inset (center): Histogram of the blue conductance trace in (a) 
computed by binning the conductance data (acquired at 40 
kHz). All histograms in (b) are obtained using a 10" 4 Go 
bin size and are normalized to the number of traces, (c-d) 
Normalized histograms of (c) step lengths (bin size: 0.004 
Go/nm) and (d) step slopes (bin size: 0.0025 Go/nm). 



bution of step-averaged conductances in traces with a 
molecular step (blue histogram in Fig. 1(b) inset) to 
that obtained by binning all data points in the traces 
(red histogram in Fig. 1(b) inset). The conductance 
peak position and half width are 0.0065 (0.0065) Go and 
± 40 (45) % for the blue (red) histograms, in excel- 
lent agreement with those for the full histogram (black). 
The slightly narrower peak width in the step-average his- 
togram indicates the peak width (red) originates primar- 
ily from variations in conductance among different junc- 
tions: Junction-elongation in individual traces evidently 
does not significantly add to the spread in conductance. 

The repeated plastic deformation of the contact region 
from measurement to measurement induces variations in 
junction structure. Understanding the impact of junc- 
tion structure on conductance is critical. We approach 



this problem in two steps. Using first-principles DFT, 
possible bonding motifs in the junction are analyzed to 
probe the selective bonding hypothesis [§]. Then, guided 
by allowed bonding motifs, 15 distinct junction geome- 
tries are examined to relate the conductance to junction 
structure. 

Our DFT calculations are performed within the gen- 
eralized gradient approximation (GGA) [17], as imple- 
mented in the SIESTA code fl8| . An optimized single-^ 
basis set is used for the Au d shell; all other orbitals are 
described by double-^ polarization basis sets [l9| . Junc- 
tions are first constructed with 6-layer Au(lll) slabs on 
either side of the BDA molecules. A supercell geome- 
try with 16 Au atoms per layer is used. All atoms in 
the junction, except for those in the bottom 3 Au layers 
of each slab, are relaxed until forces are < 0.05 eV/A. 
The conductance is obtained using a coherent scattering- 
state approach based on DFT [lj|. The BDA junction 
is divided into three regions: left bulk, center resistive 
region (including 4 Au layers on cither side of molecule), 
and right bulk. The Hartree potential in the center re- 
gion smoothly matches that of the bulk regions, which 
are infinitely repeated away from the junction to sim- 
ulate open boundary conditions. The two-dimensional 
Brillouin Zone is sampled by an 8 x 8 Monkhorst-Pack 
&)||-mesh in most junctions; in the (III', III') junctions, a 
16 x 16 mesh is found to be necessary for converged con- 
ductance. Energy- and fey -dependent scattering states 
are constructed with incoming and outgoing states deter- 
mined from the bulk Au complex band structure. The 
zero-bias conductance is computed from the Landaucr 
formula G = Go Efc,, Tr(iTi), where the transmission ma- 
trix t is evaluated at the Fermi energy, Ep- 

Unlike thiol- Au linkages [HI, we find that the local 
amine- Au bonding motif is remarkably well-defined, with 
the amine group only binding to undercoordinated Au 
atop sites. If the amine group is initially positioned at the 
bridge or three-fold site between Au atoms, the molecule 
is found to relax to the nearest atop site, or drift away 
from the Au contact. Furthermore, typical BDA bind- 
ing energies on undercoordinated Au sites (-0.4 to -0.6 
eV) are significantly larger than that on the atop site of 
atomically flat Au(lll) (-0.1 eV), in accord with previ- 
ous DFT calculations on alkanediamines Q. This pro- 
vides detailed support for the initial hypothesis 0: the 
strongly basic amine group readily donates its lone pair 
to form Au-N bonds with an undercoordinated Au atom 
where the 6s orbital is more available, in contrast to the 
relatively unreactive sites on flat Au(lll) 

We next develop model contact geometries with varia- 
tions in junction structure, considering a total of 15 fully- 
relaxed molecular junctions. These geometries differ in 
the molecular tilt angle (0°, 10°, 31°, 56° with Au sur- 
face normal), bonding configuration (cis/ trans), binding 
sites, and Au contact structures. Both symmetric and 
asymmetric junctions are examined. Six distinct Au con- 
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FIG. 2: (a) Au contact motifs, (b) side junction, (c) (III, I) 
junction. In (b), one amine group is bonded to an adatom 
(I), and the other to an atom at the base of a pyramid (III'). 
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FIG. 3: Energy-dependent transmission for the (III, I) junc- 
tion. Isosurface plots illustrate the molecular orbitals respon- 
sible for the transmission peaks. 



tact structures with varying degrees of coordination are 
considered (Fig. 2(a)). Except for the side junction in 
Fig. 2(b), all other junctions are labeled according to its 
two Au contact structures (e.g., (Ill, I) in Fig. 2(c)). The 
calculated binding energy for the BDA varies modestly 
among the 15 junction structures (0.8-1.3 eV). 

Fig. 3 shows a typical calculated transmission curve, 
T(E), at zero bias, obtained for the (III, I) junction. 
The peaks centered at -1.2, 2.7 and 3.2 eV derive respec- 
tively from the BDA highest occupied molecular orbital 
(HOMO), lowest unoccupied molecular orbital (LUMO) 
and LUMO + 1 states. The transmission for all 15 junc- 
tions is plotted on a log-scale in Fig. 4(a), and illustrates 
the remarkably narrow spread in T(E) near Ep. More- 
over, for all junctions, the transmission from ~ -1.5 to 
0.5 eV is found to be associated with a single eigenchan- 
nel arising from the BDA HOMO, consistent with 
experiment Q. For the 15 junctions considered here, the 
calculated linear response conductance GoT(Ep) has a 
mean value of 0.046 Go, and ranges from 0.032 (-31.1 
%) to 0.062 (+35.4 %) G . Notably, the spread is simi- 
lar to the half-width observed in the histogram of mea- 
sured step-average conductances (Inset, Fig. 1(b)); the 
calculated conductance magnitude, however, exceeds the 
measured value by about a factor of seven. 

The geometric feature with the greatest impact on 
T(E) is n, the total number of Au nearest neighbors of 
the two binding sites. The right tail of the resonance 
peak, centered at ~ -1 eV, is well fit to a Lorentzian 
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FIG. 4: Calculated transmission on (a) a log scale for all 
15 junctions (black: (III', III') junctions), and (b) a linear 
scale for 4 representative geometries, vertically offset for clar- 
ity. Inset: Scattering state PDOS on Au binding sites, in the 
absence of BDA. (c) Calculated conductance versus conduc- 
tance predicted from Lorentzian hts (red dashed lines in (b)). 
The resonances centered at ~ -1 eV have a maximum trans- 
mission of 0.97-0.99 Go- Numbers refer to n described in the 
text. 



form T(E) = A/((E - E ) 2 + (T/2) 2 ), for each junction 
with n > 6 (Fig. 4(b)). For n = 6, additional structure 
in T(E) appears near Ep. The relatively good fits for 
n > 6 reflect a modest electronic coupling between the 
BDA HOMO and an approximately constant density of 
states (DOS) in the leads. The structure in T(E) near 
Ep for n = 6 is due to additional features in the projected 
DOS on the Au contacts (Fig. 4(b) inset), which arise 
from the interaction of bulk and surface states with the 
s-like orbitals of the undercoordinated adatom. Overall, 
the transmission resonance energy E$ (-1.27 to -0.94 eV) 
and peak width T (0.34 to 0.56 eV) exhibit modest vari- 
ation. Further, the small increase in peak width Y with 
coordination n leads to an increase in conductance, evi- 
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dent from the comparison of the Lorentzian conductance 
(predicted from the fit) with the calculated conductance 
in Fig. 4(c). However, the extra structure in T(E) for 
ri = 6 increases the conductance beyond the Lorentzian 
model and reduces the variation of conductance with n. 

From the experimental and theoretical results, we can 
conclude that the local amine- Au bonding motif does not 
vary significantly across different junction geometries. 
Large changes in the molecular tilt angle and bonding 
configuration can be accommodated with relatively small 
variations in the Au-N-C bond angles (which range from 
112° to 130°) and N-Au bond lengths (which range from 
2.38 to 2.55 A). Because of the intrinsic isotropy of the 
Au 6s orbital, different amine- Au bond orientations have 
minimal impact on the electronic coupling. We note that 
junction-elongation process does not add significantly to 
the spread in conductance (Fig. 1). Most of the mea- 
sured conductance steps are relatively short, consistent 
with breaking the N-Au bond 21]. The well-defined 
amine- Au bonding motif, taken together with electronic 
flexibility of the amine- Au bond, result in a narrow con- 
ductance distribution across distinct junctions (Fig. 1). 

Having developed a consistent picture for junction 
structure in agreement with experiment, we may pro- 
vide a direct assessment of the DFT framework for com- 
puting conductance. Our calculations overestimate the 
conductance by ^7x. Given the relatively small effects 
of junction geometry on conductance, this discrepancy 
is unlikely to arise from finite temperature fluctuations. 
Instead, the deviation likely originates from the use of 
a mean-field potential in DFT. Indeed, many-electron 
interactions are essential to accurately position frontier 
molecular resonance energies in a junction, and the DFT 
errors can be large 22j. Direct treatment of such in- 



teractions in transport through molecules is challenging 
and to date, studies have been confined to model sys- 
tems 23, 13, 25, 2(1 27, 28]. However, many-electron in- 
teractions are expected to reduce the DFT conductance, 
shifting the molecular resonances away from Ep 2(| 27], 
and narrowing the transmission peak widths in the weak 
coupling limit 25j, |27| • To have a sense for the magnitude 
of the changes required, either reduction of the width of 
the resonance T from ~0.5 to 0.2 eV or shift of the po- 
sition of the resonance (Eo) from -1 eV to -3 eV relative 
to Ep, would bring the calculated conductance into line 
with the measured value. On physical grounds, the res- 
onance position Eo is probably too shallow because of 
self-energy errors inherent to the use of LDA or GGA 
Kohn-Sham eigenvalues as quasiparticle energies [13] . A 
more complete first-principles treatment of electron cor- 
relation effects [29} suggests that the resonance should be 
deeper [22| . Our simple estimate, including electrode po- 
larization in a physically-motivated image charge model 
[22] ] ,suggests that a resonance position of -3 eV is plausi- 
ble [30( . Since the low bias conductance is determined by 
transmission through the tail of the HOMO-derived res- 



onance, these shifts significantly alter the conductance, 
but they minimally affect the other calculated trends. 

Portions of this work were performed at U.S. De- 
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AC0298CH10886). This work was also supported by the 
Nanoscale Science and Engineering Initiative of the NSF 
(Award Numbers CHE-0117752 and CHE-0641532), the 
New York State Office of Science, Technology, and Aca- 
demic Research and by NSF Grant No. DMR04-39768. 
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